Astronomy & Astrophysics manuscript no. gradz l3nov201 1 


©ESO 2011 


November 17, 2011 





MASSIV: Mass Assembly Survey with SINFONI in VVDS^ 
III. Evidence for positive metallicity gradients in z ~ 1.2 star-forming galaxies 

J. Queyrel''^, T. Contini^-^, M. Kissler-Patig^ B. Epinat''^, P. Amram'^, B. Garilli^ O. Le Fevre"*, J. Moultaka^-^, L. 
Paioro^, L. Tasca^, L. Tresse'^, D. Vergani^, C. Lopez-Sanjuan^, and E. Perez-Montero' 

' Institut de Recherche en Astrophysique et Planetologie (IRAP), CNRS, 14, avenue Edouard Belin, F-31400 Toulouse, France 

^ IRAP, Universite de Toulouse, UPS-OMP, Toulouse, France 

^ ESO, Karl-Schwarzschild-Str.2, D-85748 Gai'ching b. Munchen, Germany 

* Laboratoire d' Astrophysique de Marseille, Universite de Provence, CNRS, 38 rue Frederic Joliot-Curie, F-13388 Marseille Cedex 
13, France 

5 lASF-INAF, Via Bassini 15, 1-20133, Milano, Italy 

* INAF-Osservatorio Astronomico di Bologna, via Ranzani 1, 1-40127, Bologna, Italy 
' Institute de Astrofi'sica de Andaluci'a - CSIC Apdo. 3004 E- 18080, Granada, Spain 

Received ; accepted 

ABSTRACT 

Context. A key open issue for galaxy evolution and formation models is the understanding of the different mechanisms of galaxy 
assembly at various cosmic epochs. Spatially-resolved spectroscopy of high-redshift galaxies has proven to be a powerful observing 
technique to address this issue. 

Aims. We aim at deriving the global and spatially-resolved metal content in high-redshift galaxies. This allows to study, mainly via 
the determination of radial abundance gradients in the disks, their star formation history and their interaction with the surrounding 
intergalactic medium. 

Metliods. Using VLT/SINFONI integral-field spectroscopy of a first sample of 50 galaxies at z ~ 1.2 in the MASSIV survey, we are 
able to measure the Ha and [Nn]6584 emission lines. Using the N2 ratio as a proxy for oxygen abundance in the interstellar medium, 
we measure the metallicity of the sample galaxies. We develop a tool to extract spectra in annular regions of these galaxies, leading 
to a spatially-resolved estimate of the oxygen abundance in each galaxy. 

Results. We derive a metallicity gradient for 26 galaxies in our sample and discover a significant fraction of galaxies with a "positive" 
gradient. Using a simple chemical evolution model, we derive infall rates of pristine gas onto the disks. 

Conclusions. Seven galaxies display a positive gradient at a high confidence level. Four out of these are interacting and one is a chain 
galaxy. We suggest that interactions might be responsible for shallowing and even inverting the abundance gradient. We also identify 
two interesting correlations in our sample: a) galaxies with higher gas velocity dispersion have shallower/positive gradients; and b) 
metal-poor galaxies tend to show a positive gradient whereas metal-rich ones tend to show a negative one. This last observation can 
be explained by the infall of metal-poor gas into the center of the disks. We address the question of the origin of this infall under the 
influence of gas flows triggered by interactions and/or cold gas accretion. 
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1. Introduction 

An key open issue for galaxy evolution and formation models 
is the understanding of the different mechanisms of galaxy as- 
sembly at various cosmic epochs. In this context, the gas-phase 
and stellar metallicities have proven to be important parameters 
to constrain the star formation history of galaxies. 

A relation between galaxy mass and metallicity, first dis- 
covered by Lequeux et al. ( 1979^ in irregular galaxies, exists 
for star-forming galaxies both in the nearby universe ( |Tremonti| 
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* This work is based on observations collected at the European 
Southern Observatory (ESO) Very Large Telescope, Paranal, Chile, as 
p^F of the Pr ograms 1 79.A-0823, 78.A-0177, and 75.A-0318. This 
work also benefits from data products produced at TERAPIX and 
the Canadian Astronomy Data Centre as part of the Canada-France-, 
HawaiiTelescope Legacy Survey, a collaborative project of NRC and 



CNRS. 



^TaLjllOO ?, 'Lama reille eTaLjllOO?]!, a nd at high redshifts 



0.7 - 3 (Lamareill e et al.| 12006 712009; Perez-Monter o et aL| 
2009; C^ueyrel et al. 2009; Erb et al. 2006; Mai ohno et al.|2008| 



Mannucci et al._20 09). This trend for more massive galaxies to 
have a higher gas-phase metallicity can be explained by various 
related properties: gas inflows and/or outflows, strength of the 
gravitational potential, efficiency of the star formation depend- 
ing on stellar mass, etc. Their relative importance is currently 
being explored in numerical simulations. 

Studying the abundance in Hii regions of nearby spi- 
ral galaxies has unveiled metallicity gradients in the local 
Universe (|Pagel & Edmunds||1981| |Vila-Costas & Edmunds 



T992[ [Considere et al."2000; Pilyugi n et al.||2004| l. hi z ~ 



galaxies, the metallicity generally decreases from the center to 

the outskirts (up to A^Z 0.1 dex kpc"'). The strength of these 

gradients seems to correlate with various parameters such as ro- 
tation velocity, morphological type, luminosity, and the presence 
of a bar. Different physical processes have been proposed to ex- 
plain these patterns, including radial gas-flows ( jKoeppen|1994| l, 
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self-regulati ng star formation ([Phillipps & Edmunds|199T] l, vari- 
able yields ("Vila-Costas & Edmunds'1992'), or continuous infall 
of gas onto the disk ( Magrini et al. 2007) . 

Merging events - which are believed to play a substantial 
role in galaxy evolution (eg. de Ravel et al. 2009 Lopez-Sanjuan 
et al.|20lT[) - seem to be a key physical process in shaping the 



metallicity gradients of interacting galaxies ( Rupke et al.|2010 i. 
Recent observations have suggested that galaxies involved in 
merging events show lower nuclear metallicities due to the infall 
of pristine gas into the nucleus. Merging events could account 
also for outliers to the mass-metallicity relation ( Kewley et al. 
2006, Michel-Dansac et al. 2008; Peeples et al. 2009, Queyrel 
et al.,|2009, ,Alonso et al.||2010; ,Montuori et al.n2010, ,Kewley 



et al.|2010| l. 



Metallicity gradients have been observed in the stellar pop- 
ulations of early-type galaxies as well. They appear to correlate 
with macroscopic properties such as, for example, stellar mass 
(eg. |Spolaor et al.| 2009'). It has been demonstrated that a metal- 
licity gradient in the stellar populations could survive a major 
merger event, although it will be weakened. These observations 
contrast the predictions of the simple monolithical collapse sce- 



nario ([Wh ite 1980] |Bekki & Shioya|1999j|Kobayashi|2004[ p 
[Matteo et al.|2009| . 

In nearby galaxies, metallicity gradients can be inferred 
from the observation of different emission-line ratios (in H ii re- 
gions), or from stellar absorption lines. In the distant universe 
(z > 0.5), low signal-to-noise ratio (hereafter SNR) and low 
spatial-resolution data make these measurements much more 
challenging. However, thanks to powerful integral-field spec- 
trographs mounted on the largest telescopes, it is now possible 
to determine spatially-resolved physical parameters of high red- 
shift galaxies. 

Many studies have taken advantage of this technique to study 
the gas dynamics of distant galaxies, from z ~ 0.6 ([Flores et al" 
2006) up to z > 2 d Forster Schreiber et al.||2009l [L aw eTaL 
2009; Epinat et al. 2009; Gnerucci et al. 2011). Recently, the 
Mass Assembly Su rvey with SINFONI in VVDS (MASSIV; 
Contini et al. 201 l| l collected data of z ~ 1-2 galaxies ob- 



served with the integral-field spectrograph SINFONI at the VLT 
( Eisenhauer et al.||2003[ l. Depending on the galaxy redshift, the 
bright Hot emission line is targeted in the J or H band. The nitro- 
gen line [Nn]6584 close to Ha allows to estimate the gas-phase 
metallicity via the N2 ratio (N2 = log ([Nn]6584/Ha)) and the 
corresponding abundance calibration ( Kewley & Dopital|2002 



[Denicolo et"al.,2002; Perez-Montero & Contini|2009| ). 

The aim of this paper is to investigate - for the first time 
at a redshift around 1 - metallicity gradients of 50 star-forming 
galaxies in the MASSIV sample. In order to achieve this goal, we 
develop a dedicated program to analyse the low SNR data cubes. 
Two companion papers ( Vergani et al.|201 1] Epinat et al.|20I l| l 
discuss the associated fundamental scaUng relations (e.g. the 
Tully-Fischer relation for disks) and the kinematical proper- 
ties of the sample galaxies. The results of these three studies 
are cross-correlated in order to investigate possible relations be- 
tween the metallicity gradients and global properties of galaxies 
atz ~ 1.2. 

The paper is organized as follows. In Section |2] we briefly 
summarise the properties of the galaxy sample, the data re- 
duction technique, the emission-line measurements, the level of 
Active Galactic Nucleus (AGN) contamination, the kinematics, 
and the mass (both stellar and dark matter halo) estimates. In 
Section|3]we present and discuss the determined metallicity gra- 
dients. Our conclusions are drawn in Section|4] 



Throughout the paper, we assume a ACDM cosmology with 
Q„ = 0.3, Qa = 0.7 and Hq = 70 kms^^Mpc"'. 



2. Data Description 

2.1. Target selection, observations and data reduction 

The galaxy sample studied in this paper is the "first epoch" 
sample of the MASSIV project (ESO Large Program , PI.: T. 
Contini). A full description of the sample can be found in Contini 
et al. ( 201 1} . We briefly summarise some properties of this sam- 
ple of 50 galaxies below. 

The galaxies were selected from the VIMOS VLT Deep 
Survey (VVDS) in the deep (Iab < 24; |Le Fevre et"ar][2005) , 
ultra-deep {Iab < 24.75; Le Fevre et al., in prep ) RA = 02h, and 
wide {Iab < 22.5; |Garilli et al.|2008| l RA = 14h,22h fields. The 
galaxies were chosen according to their redshift such that their 
Ha line (or [Oni]5007 in a few cases) was visible in the J or 
H band, and was not affected by a bright OH sky-line. Galaxies 
were selected to be star-forming on the basis of their [On] 3727 
emission line strength. The observations have been performed 
between April 2007 and August 2009. The general properties 
(RA, DEC, redshift, /-band magnitude) and characteristic ob- 
servational configuration (spectral band, adaptive optics or not, 
exposure time, spatial resolution) of each galaxy are gathered in 
Table [T| Most (85%) of the galaxies in this "first epoch" sample 
have been observed in a seeing-limited mode (with a spatial sam- 
pling of 0.125")- However, seven galaxies have been acquired 
with Adaptive Optics assisted with a Laser Guide Star (AO/LGS, 
0.05" spatial sampling). Among the 50 "first epoch" MASSIV 
galaxies observed with SINFONI, four galaxies have not been 
detected (see |Epinat et al. 201 1 , for details). 

The data reduction was performed with the ESO SINFONI 
pipeline, using the standard master calibration files provided by 
ESO. The absolute astrometry for the SINFONI data cubes was 
derived from nearby bright stars also used for PSF measure- 
ments. Custom IDL and Python scripts have been used to flux 
calibrate, align, and combine all the individual exposures. For 
each galaxy a non sky-subtracted cube was also created, mainly 
to estimate the effective spectral resolution. For more details on 
data reduction, we refer to |Epinat et al.| ( |20TT] l. 

/-band images for all the galaxies were obtained through 
CFHT Megacam imaging (for the 22h and 02h field, from the 
CFHTLS "best seeing") and CFHT-12K imaging (for the 14h 
field). These images were used for two purposes: (i) refinement 
of the SINFONI astrometry, using the relative position of the 
PSF star, and (ii), deriving morphological parameters (used as 
inputs to the kinematics modeling, see , Epinat et al.,201 1) . 



2.2. Measurement of emission lines 

As mentioned in the previous section, galaxies were selected to 
be star-forming such that emission lines from Hii regions are 
visible in the near infrared. Given the redshift range of MASSIV 
galaxies, we targeted, in most of the cases, the rest-frame optical 
Ha emission line in the J (for z < 1 .2) or // (for z > 1.2) band. 
For four objects (VVDS020 126402, 020167131, 020306817, 
and 220148046) we observed the [Oiii]5007 fine which allowed 
to increase the redshift window (two of them were not detected). 
These latter galaxies are not considered in this paper as we 
cannot derive any metallicity from the [Om]5007 emission-line 
alone. In the former cases, we took advantage of the proximity 
of the nitrogen line to Ho- to derive the metallicity of a majority 
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Table 1. General properties and SINFONI observing log of the 50 "first epoch" galaxies of MASSIV. Column 1: identifier in the 
VVDS catalog. Columns 2 & 3: RA & DEC in degrees. Column 4: spectral band of the observation. Column 5: seeing-limited 
mode (0) or Adaptive Optics (1), Column 6: on-source total exposure time. Column 7: /-band magnitude in the AB system. 
Column 8: Mean spatial resolution (PSF) as measured on a nearby star. Column 9: redshift from Ha (or [Oin]5007 when marked 
with an asterisk) as measured with SINFONI. For the four galaxies (VVDS020126402, VVDS020217890, VVDS020306817, and 
VVDS220071601) which have not been detected with SINFONI, redshift comes from VIMOS spectra 
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of the galaxies using the [N n]6584/HQ' line ratio as a proxy for 
oxygen abundance (see sect.|3]l. 



2.2.1. Integrated fluxes 

The global integrated flux of each emission line was calculated 
using a mask designed for each galaxy. The flux was collected 
over a region where the signal-to-noise ratio of the Ha line is 
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above a fixed threshold of SNR > 2. The SNR maps were pro- 
duced using a Gaussian-fitting procedure written in IDL, already 
used and described in Epinat et al. ( 201 1| ). A Gaussian spatial 
smoothing of 2 x 2 spaxels (0.25"x 0.25") was used, increasing 
the SNR without degrading the final spatial resolution (~ 0.65" 
on average). The line position map produced by this procedure 
was used to shift in wavelength every single spectrum such as 
to compensate for the Doppler effect corresponding to the mea- 
sured Ha velocity at that position. In this way, the broadening of 
the line by large-scale rotation in the spatially-integrated spec- 
trum was largely removed and the line profile more easily fitted 
with a Gaussian, increasing the accuracy of line measurements. 

The Ha fluxes were measured on flux-calibrated integrated 
spectra using a Gaussian fit and a flat continuum. These fluxes 
have been corrected for dust reddening using the extinction co 
efficient derived from the SED fitting (see Contini et al.||2011 



Table 2. Ha flux, A^2 - log([Nii]6584/Ha) emission-line ratio, 
and Ha-based star formation rates (corrected - and not - for dust 
reddening) for each galaxies in our sample. 



|Vergani et al.||20lT| l. The line ratios, on the other hand, were 
determined from the integrated spectra in counts. Indeed, the 
flux calibration procedure would have added noise on top of 
the akeady low SNR data. Further, the lines of interest (namely. 
Ha and [Nn]6584) are close enough in wavelength to assume 
that the sensitivity curve is constant over the wavelength range 
of interest. The same argument justifies the fact that we did 
not correct the emission-line ratios for differential extinction. 
Table |2] lists our final integrated Ha flux, N2 emission line ra- 
tio (A^2 = log([Nii]6584/Ha) and star formation rates (SFR, 
corrected or not for dust reddening). Among the 44 galaxies de- 
tected in Ha with SINFONI, we have been able to measure the 
[Nii]6584 emission-line in 34 galaxies and hence derive inte- 
grated N2 emission-line ratio and metallicity for these objects. 

2.2.2. Spatially-resolved flux and line ratio 

For many galaxies in our sample, the low SNR of the data cubes 
prevents us from measuring in each spaxel, with high confi- 
dence, the emission lines surrounding Ha : namely [Nn]6584, 
[S ii]6717,31, and [Oi]6300 which is detected only in one 
galaxy: VVDS 140096645. These lines might either not be bright 
enough to be detected, or a bright OH sky-line is too close. To 
increase the SNR of our final spectra, we summed the spectra of 
the spaxels gathered in a specified region, and measured the line 
fluxes in the resulting spectra. 

We developed a specific Python procedure to do this analy- 
sis within our data cubes. The program allowed to define a re- 
gion spaxel-by-spaxel, or to define a region with contours on 2D 
maps (such as Ha flux, SNR, etc maps). In a given region, the 
spaxel-to-spaxel line shift due to the rotation velocity of the gas 



along the line-of-sight was fully corrected for (see §2.2.1 1. The 
integrated spectrum in a defined region was then fitted with a flat 
continuum and two Gaussians, one for Ha and one for [N ii]6584 
(the sulfur doublet was neglected for this analysis as it was de- 
tected in a few cases only). The width of each Gaussian was set 
to be the same as the emission of the collisional and recombi- 
nation lines traces the same ionised gas in the galaxies. The fit 
was weighted using the corresponding non sky-subtracted spec- 
trum, i.e. giving a lower weight to the channels affected by the 
subtraction of a strong sky-line. We then estimated the Icr er- 
ror on the fluxes with a Monte-Carlo technique, in which the fit 
was considered to be a "noise-free" model, to which we added 
a Gaussian noise with a standard deviation corresponding to the 
residuals from the initial spectrum minus the "noise-free" model. 
This operation was repeated a hundred of times, which gave a set 
of parameters on which the Icr deviation was computed. Fig.[T] 
shows a result of our fits in different regions for two galaxies. 



Galaxy 


F(Hq') 




N2 


S FRlia 






[10""ergs-'cm--] 






[Mo yr'] 


[Mo yr'] 


(1) 


(2) 




(3) 


(4) 


(5) 


020106882 


14.0 ±0.8 


-0.71 


±0.12 


13.3 


38.1 


020116027 


12.4 ±2.9 


-1.02 


±0.15 


14.7 


42.2 


020126402 












020147106 


46.8 ±5.8 


-0.95 


±0.14 


54.3 


92.1 


020149061 


20.6 ± 1.6 






16.0 


45.9 


020164388 


24.1 ± 1.3 


-0.86 


±0.10 


21.1 


46.5 


020167131 












020182331 


15.2 ±5.0 






10.5 


51.1 


020193070 


6.7 ± 0.9 


-0.56 


±0.13 


3.0 


14.4 


020208482 


2.6 ±0.5 






1.2 


2.6 


020214655 


23.3 ± 1.1 


-0.68 


±0.06 


10.6 


51.7 


020217890 












020239133 


7.8 ± 1.0 






3.4 


16.5 


020240675 


8.5 ± 1.0 






7.1 


15.7 


020255799 


5.8 ±0.7 


-0.65 


±0.15 


2.6 


12.8 


020261328 


9.5 ± 0.9 


-1.06 


±0.48 


11.2 


19.0 


020278667 


2.4 ± 1.3 






1.1 


4.2 


020283083 


13.4 ± 1.0 


-0.70 


±0.12 


10.3 


17.4 


020283830 


11.0± 1.1 


-0.82 


±0.17 


10.3 


50.4 


020294045 


15.5 ±2.8 






6.5 


14.3 


020306817 












020363717 


37.4 ±4.3 


-0.87 


±0.09 


31.5 


53.4 


020370467 


17.3 ±2.9 


-0.69 


±0.11 


14.6 


71.2 


020386743 


17.4 ± 1.7 


-0.91 


±0.08 


8.1 


39.4 


020461235 


7.5 ±2.0 


-0.67 


±0.12 


3.3 


9.6 


020461893 


14.5 ± 1.2 






6.7 


14.9 


020465775 


14.5 ± 2.8 






12.8 


62.2 


140083410 


47.0 ±3.3 


-0.72 


±0.08 


16.8 


37.1 


140096645 


55.4 ±4.2 


-0.25 


±0.02 


20.9 


102.1 


140123568 


5.8 ± 1.2 


-0.69 


±0.10 


2.4 


11.8 


140137235 


9.6 ±2.3 


-0.79 


±0.13 


4.4 


21.6 


140217425 


104.7 ± 4.9 


-0.44 


±0.02 


41.0 


200.0 


140258511 


42.5 ± 2.7 


-0.51 


±0.10 


30.0 


146.4 


140262766 


13.1 ±5.0 


-0.75 


±0.15 


10.0 


10.0 


140545062 


29.4 ±3.3 


-0.75 


±0.06 


13.4 


22.7 


220014252 


50.3 ±3.1 


-0.67 


±0.08 


40.5 


197.5 


220015726 


47.5 ± 6.7 


-0.69 


±0.04 


37.0 


106.4 


220071601 












220148046 












220376206 


72.4 ±5.1 


-1.04 


±0.10 


51.2 


249.7 


220386469 


17.1 ±2.4 


-1.14 


±0.16 


7.5 


36.5 


220397579 


65.1 ±9.1 


-1.07 


±0.07 


29.4 


143.2 


220544103 


55.0 ± 1.9 


-0.56 


±0.09 


52.2 


117.5 


220544394 


24.3 ± 1.7 


-0.89 


±0.07 


10.3 


50.1 


220576226 


31.2 ± 1.6 


-0.69 


±0.04 


13.6 


66.3 


220578040 


20.5 ± 1.9 


-0.80 


±0.07 


9.4 


20.8 


220584167 


64.3 ± 1.3 


-0.91 


±0.10 


68.6 


202.6 


220596913 


30.0 ± 1.0 






23.2 


36.1 


910193711 


32.9 ± 18.4 


-0.68 


±0.09 


40.6 


197.7 


910279515 


14.9 ±2.1 


-0.83 


±0.21 


14.2 


69.0 



Finally, we estimate spatially-resolved Ha and [Nii]6584 
emission-line flux for 26/34 (~ 75%) galaxies in the sample. The 
remaining 8 galaxies have either a too low SNR or the [N ii]6584 
emission-fine is polluted by sky-fine residuals. 



2.3. AGN contamination 



A recent study (Wright et al. 2010 1 has shown that unveiling 
the presence of AGN in high-redshift galaxies is a difficult ex- 



J. Queyrel et al.: MASSIV: Mass Assembly Survey with SINFONI in VVDS 



5 





' Region 1 










40 

t "A' *" 









Fig. 1. Examples of spatially-resolved emission-lines measure- 
ments in two MASSIV galaxies (Top: VVDS220578040, bot- 
tom: VVDS140217425). Color-coded images show the regions 
defined with Ho- contours used to derive the metallicity gradi- 
ents. The ID spectra integrated over the different regions are 
shown with solid black lines. Model (in red) and sky (dashed 
line) spectra are also shown. The sky spectrum is used as a 
weight in the fitting procedure. Region 1 is the outer most with 
the region number increasing towards the center. 



ercise. In the case of metallicity studies, in which abundances 
are deduced from the ratio of diff'erent emission lines of the 
ionised gas, it is critical to check that the intensity and width 
of these lines are due to star formation and not related to any 
non-thermal nuclear activity. The common way to disentangle 
AGN contribution from star-forming galaxies consists in com- 
paring the relative intensity of the main nebular emission lines 
(mainly [Oiii]5007, Hy8, Ha, and [Nii]6584) in a diagnos tic dia- 
gram, so-called BPT diagram ( [Baldwin et al.|I981 ; Kewley et al.| 
200 1 [ l. Various physical conditions in the ISM - SFR, ionisa- 
tion parameter, metallicity and/or chemical composition - have 



been invoked to explain the fact that some high-z star-forming 
galaxies lie in the transition region of the local BPT diagram 
(as defined by the SDSS galaxies) between star-forming galax- 
ies and AGN hosts. Wright et al. (2010 1 have been able - thanks 
to high resolution adaptive optics observations - to subtract the 
active nuclear emission in a z ~ 1.6 galaxy (HDF-BMZ1299), 
and have shown that the residual extended star-forming emission 
was characteristic of a local SDSS star-forming galaxy, whereas 
the integrated emission would have placed the object in the tran- 
sition region. 

The nature of our observations did not give us simultane- 
ous access to the set of emission lines ([Oiii]5007, Hfi, Ha, and 
[N ii]6584 or [S n]6717,31) commonly used in standard diagnos- 
tic diagrams. However, for all but two objects in our sample, the 
emission-line ratio N2 = log([N n]6584/HQ') is lower than -0.5, 
with a median N2 value of -0.72. Such low values are indica- 
tive for a very low contamination by AGN in our sample (eg. 
[Baldwin et al.|1981| l. 




-1 -0.5 



Fig. 2. Left: Distribution of the N2 ratio in the MASSIV "first 
epoch" sample. The black line represents the ratio integrated 
over the whole galaxy. The red line is the N2 ratio measured 
in the nuclear region. Right: N2 global ratio versus N2 nuclear 
ratio. 



For 24 galaxies of our sample we calculated, following 
[Wright et al.j ( [2010| l, the N2 "concentrated ratio", correspond- 
ing to the value in the nuclear region of the galaxy. We defined 
the nuclear region as the spaxel with the highest Ho- flux along 
with its 8 nearest neighbours (corresponding to a 0.7" diame- 
ter aperture, matching our mean spatial resolution of 0.65"). We 
assumed this aperture to be small enough to probe the inner nu- 
cleus part as objects usually span from 1" to 2" with SNR > 2 
in our observations. 

Fig. |2] shows i) the distribution of the 24 galaxies as a func- 
tion of their global and nuclear N2 ratio (left panel) and ii) 
the relation between the global and nuclear N2 ratios for each 
galaxy (right panel). The median values of N2 for each distri- 
bution are not very different (A ~ -0.07). The median nuclear 
N2 ratio is lower than the global ratio which would not be the 
case if a significant fraction of our sample galaxies were host- 
ing an AGN. When comparing the global ratio to the nucleus 
ratio distribution, the highest bin does not shift and contains a 
single object . We investigated in more detail the galaxy in this 
bin: VVDS 140096645. It shows the following high N2 ratios: 
^^global ~ -0.252 and A^2jjm.jg^ - -0.294. Looking further 
into its integrated spectrum (global and nuclear, see Fig. [3]), we 
noticed that: (i) the emission lines are broad, which is a possible 
sign of nuclear activity, (ii) the two nitrogen lines are clearly vis- 
ible, as is the sulfur doublet, and the [Oi]6300 line, which alto- 
gether are characteristics of LINER galaxies, often associated to 
AGNs ( [Heckm an 1980) - or violent episodes of star-formation 
in high metalUcity galaxies ( [Terlevich & Melnick|1985 1. 

The galaxy with the next highest global N2 value is 
VVDS 1402585 11 (A^2gjob^j = -0.509 and A^2jj^^.jg^ = 
-0.345). Its global and nuclear spectra show no obvious features 
of nuclear activity, and the N2 ratio are both < -0.3. 

We conclude that our sample does not suffer from significant 
AGN contamination. The only candidate for nuclear activity is 
VVDS 140096645. Similar conclusions were reached based on 
composite ID VIMOS and SIN FONI spectra of the MASSIV 
sample (see Contini et al.|201 1 1. 



2.4. Kinematics and close environment classification 

A kinematics and close environment classification of the 
MASSIV galaxies has been performed and is described in de- 
tails in [Epinat et al.[ ( [2011| l. This classification is based on the 
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Fig. 3. Integrated spectra of VVDS 140096645. Top: spectrum 
of the disk (total integrated with central contribution removed). 
Middle: spectrum of the central region surrounding the peak of 
Ha flux. Bottom: spectrum from the whole spatial extent. From 
left to right, the following emission lines appear: [Oi]6300, 
[Nii]6548, Hot, [Nii]6584, [Sii]6717, [Sn]6731. The red and 
blue lines are least-square fits to the data. 



/-band morphology, the Ha flux maps, the Ha velocity fields 
and their modeling assuming a rotating disk. 

In the present paper, we exploit two types of classes: a) the 
dynamical state of galaxies (rotating or non-rotating), and b) the 
close environment of galaxies. The first class relies i) on the 
agreement between the position angles of the major axis deduced 
from the morphology and from the kinematics, and ii) on the ac- 
curacy of the rotating disk model. In the case the average SNR is 
lower than 5, this classification is believed to be not reliable and 
is thus not used. The second class relies on the detection of com- 
panions in both /-band image and Ha maps at a similar redshift 
than the main source. In some cases, the kinematics maps (ve- 
locity fields and velocity dispersion maps) suggest the presence 
of some companions about merging or along the line-of-sight. A 
confidence flag was assigned to the environment classification, 
ranging from A (/,90%, confident) to C (~ 50%). 

We used these classifications to help us constraining the ori- 
gin of the metallicity gradients in our MASSIV star-forming 
galaxies. 

2.5. Stellar and dark matter halo masses 

The stellar masses used in our study were obtained using 
the SED fitting technique. Stellar population synthesis mod- 
els ( |Bruzual & Ch ariot 2003) were used to match photomet- 
ric and spectroscopic data from our MASSIV sample, using the 
GOSSIP tool ( [Franzetti et al.|2008[ ). A |Salpeter| ( [T955] l IMF was 
assumed, stellar formation timescales and the extinction param- 
eter E{B - V) were allowed to range from 0. 1 to 1 5 Gyr, and from 
to 0.3, respectively. The GOSSIP tool returned the best-fit pa- 
rameters as well as a Probability Distribution Function (PDF) for 



Fig. 4. Comparison between diff'erent calibrations used to derive 
the oxygen abundance from the N2 parameter All these calibra- 
tions are consistent within their intrinsic dispersion (~ 0.3 - 0.4 
dex), especially in the range of N2 parameter (-1.0 < N2 < 
-0.5) covered by MASSIV galaxies (see Table|2]and Figure|2]). 



each of them, following Walcher et al.|(|2008j). The median and 
the standard deviation of the PDF were used to recover the pa- 
rameter estimates and their associated errors. These latter values 
were used for this study. The procedure is described in detail in 
IContini et aL] ( |20TT] l. 

The dark matter (DM) halo masses (Mhaio) used hereafter 
were computed using a spherical virialized collapse model 
dPeeblesI 19801 [White & Frenk|1991i|Mo & White|2002D : 



(1) 



where § is the universal gravitational constant, z is the red- 
shift of the galaxy and Vmax is the maximum rotational velocity 
computed in Epinat et al. (2011 1. In order to compute associ- 
ated uncertainties, the Monte Carlo method used in'Epinat et al. 
(201 1) for the uncertainties on Vmax was extended until the halo 
mass. The contributions of the two sources of uncertainty on the 
velocity (inclination and modeling) were added quadratically to 
compute the final uncertainty. 

Equation [T] makes the assumption that the plateau has been 
reached and that y,nax traces the halo circular velocity. However, 
this assumption is probably not correct for non-rotating systems. 
The stellar and DM halo masses derived for our galaxies are 
listed in Table |3] 



3. Metallicity gradients 

In the local Universe metallicity gradients in spiral galaxies are 



commonly meas ured ([Pagel & Edmunds 



Edmunds] [1992} [Considere et al.||2000[ [Pilyugin et al.||2004| l 



1981 Vila-Costas & 



These gradients are generally negative (metallicity decreasing 
from the center to the outer parts) and their amplitude is typically 
> -0.1 dex kpc ' as traced by the metallicity in Hii regions. In 
the Milky Way, a gradient of -0.07 dex kpc"' is observed. 
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Several physical processes can be responsible for such gradi- 
ents ( Goetz & Koeppen|1992 l. Radial gas flows draining metal- 
rich interstellar gas from the outer parts of the galaxy into the 
center are believed to play a key role ( [Tinsley & Larson|1978| l. 
It requires that the infall timescale of gas onto the disk is faster 
than the star formation timescale. Indeed, [Koeppen ( 19941 has 
shown that the presence of both infall of pristine gas and ra- 
dial flows into the disk is very efficient in creating an abun- 
dance gradient. Phillipps & Edmunds ( 1991 1 further suggested 
that self-regulating star formation rates varying with the galac- 
tocentric distance could generate such gradients. In contrast, for 
the case of interacting galaxies, Rupke et al.' (2010) claim that 
radial flows of low-metallicity gas from the merging galaxies 
can explain the low oxygen abundances observed at their cen- 
ter In these cases, the radial mixing of gas may flatten existing 
metallicity gradients. 

At high redshifts, metallicity gradients are harder to detect as 
the gas-phase metallicity of individual galaxies can only be mea- 
sured with a limited accuracy. Collisional emission lines from 
the metals and recombination lines from hydrogen and helium 
are the only indicators allowing to estimate the oxygen abun- 
dance in the gas-phase, needed for a direct comparison with the 
local galaxies. 

In the present study, we used the ratio of the [Nn]6584 ni- 
trogen line to the Ha B aimer line as a proxy for the oxygen 
abundance. Computing this particular ratio has two practical ad- 
vantages: the [Nii]6584 and Ha emission lines are very close in 
wavelength, so that i) the differential extinction due to dust at- 
tenuation can be neglected; and ii) the relative flux calibration 
can be considered to be constant over the spectral range. To de- 
rive from this ratio the gas-phase oxygen abundance, we used 
the calibration proposed by Perez-Montero & Contini ( 2009[ l: 



U+log(0/H) = 9.07 +0.79 xN2, where N2 = log 



[Nn]6584 



(2) 



This calibration was computed using emission-line objects (star- 
forming galaxies and H ii regions) in the nearby universe with 
an accurate oxygen abundance obtained from the measurement 
of the electronic temperature. The calibration has an intrinsic 
scatter of 0.34 dex, mainly due to the second-order dependence 
of the N2 parameter on the ionisation parameter and on the 
nitrogen-to-oxygen abundance ratio. As shown in Figure |4] the 
PMC09 calibration of the N2 parameter used in this analysis is 
totally consistent, for both the low- and high-metallicity regime, 
with other calibrations of the same parameter found in the lit- 
erature. Although the N2 parameter was calibrated by PMC09 
using data with a determination of O/H abundance based only 
on the "T-method", and hence not very well defined for the 
high-metallicity regime, the PCM09 calibration agrees very well 
within its dispersion limits (about 0.3 - 0.4 dex) with other re- 
lations that are partially based on photoionization models or 
strong-line determinations of the metallicity (Denicolo et al. 
2002] [Pettini & P agel 2004; Nagao et al.1|2006 ). This is espe- 
cially true in the range of N2 parameter (-1.0 < N2 < -0.5) 
covered by MASSIV galaxies (see Table|2]and Figure|2]). 

3.1. Observed metallicity gradients 



As reported in ^2.2[ we measured Ha and [Nii]6584 line fluxes 
in spatially-resolved regions for 26 galaxies. These regions, de- 
fined by Ha isoflux contours, are centered around the Ha peak 
in each galaxy, which most often corresponds to the kinematical 
center of the galaxy. The widths of the annular regions were ad- 
justed such as to have a high enough SNR for the measurement 



of the [Nn]6584 line and, consequently, the N2 ratio. This al- 
lowed to quantify the radial behaviour of the Ha and [Nii]6584 
lines, i.e. of the gas-phase metallicity from the inner to the outer 
parts of each galaxy. 

The derived metallicity gradients and the integrated metal- 
licities are listed in Tab [3] The metallicity is estimated as a func- 
tion of radius — each region having a mean radius and a metal- 
licity (uncertainties on metallicity estimates are dominated by 
measurements errors). The radius of a region con^esponds to the 
mean radius between its outer and inner contours. The radius of 
a contour is approximated by the radius of a perfect circle with 
the same perimeter Depending on the galaxy inclination, the in- 
trinsic deprojected radius can be underestimated up to a factor 
of ~ 1.5. Metallicities of the different regions in each galaxies 
are plotted in Fig.|5] along with the best fit line to the data. The 
figure shows that the detected radial gradients are in general very 
weak, some being positives and other negatives. 

Fig. |6] (left panel) shows the distribution of metallicity gra- 
dients (in dex kpc ') for the 26 MASSIV galaxies, for which 
we were able to measure [N ii]6584 in different regions. The his- 
tograms show the full sample (in black), the isolated galaxies 
(in red) and the interacting ones (in dashed blue), as classified 
in §2.4 For the two latter distributions, we only considered ob- 



jects classified with a high or medium confidence level (A or B), 
excluding two isolated galaxies classified with a low confidence 
level C (VVDS220578040 and VVDS910193711). In the major- 
ity of the cases, no clear gradient is detected: the median value 
of the total sample is 0.0040+0.037 dex kpc"' . The distributions 
for isolated (16 objects) and interacting (8 objects) galaxies show 
similar shapes, with median values 0.0018 + 0.036 dex kpc"' for 
the isolated galaxies, and 0.020 + 0.038 dex kpc"' for interacting 
ones. 

Interestingly though, for twelve of the 26 galaxies, a gradi- 
ent is detected with a > Icr confidence and with nearly the same 
proportion of positive (seven) and negative (five) gradients (see 
Fig.|5]l. Thus, contrary to the global trend in the local Universe, 
where the gas-phase metallicity of disk galaxies generally de- 
creases with galactocentric radius, seven of our galaxies have 
larger metallicities in the outskirts than in the center (where the 
center is defined as the maximum of the Ha flux). 

About 1/4 of our sample galaxies thus show a positive metal- 
licity gradient at a > Icr confidence level (see blue area in 
Fig. |5]l. Of these, VVDS220376206 displays a positive gradi- 
ent with 2cr confidence, while in the galaxies VVDS 140217425 
and VVDS220014252 a positive gradient of ~ 0.02 dex.kpc"' 
is detected with > 5cr confidence. Among the seven galax- 
ies, four are classified as interacting systems, while the three 
other are isolated. One of the isolated galaxies, the secure case 
VVDS 1402 17425, appears as a "chained-galaxy", with large 
clumps in the outskirts that could be interpreted as minor merg- 
ers. About half (4/8) of the interacting galaxies have a positive 
metallicity gradient while it concerns only 3/16 (~ 20%) of the 
isolated ones. We thus tentatively conclude that the majority of 
the galaxies showing a positive metallicity gradient are interact- 
ing. It is interesting to note also that among the four interacting 
systems, one galaxy only is classified as a rotating disk. 

On the other hand, five galaxies display clear negative metal- 
licity gradients at a > Icr confidence level (of which one at > 2cr 
and three at > 3cr confidence level; see Fig. |5 and Table |3]l. 
Among these five galaxies, three are isolated and two are inter- 
acting. Four of these five galaxies are classified as rotating disks. 

Although rare, this is not the first time that positive gradi- 
ents of oxygen abundance were found. Recently, Werk et al. 
( |2010| l reported a positive gradient in a local galaxy and pro- 
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Fig. 5. Metallicity gradient for the 26 MASSIV galaxies with spatially-resolved metallicities. The x-axis represents the mean radius 
(in kpc) of the region relative to the Ha center. The y-axis indicates the corresponding metallicity in dex (12 + log{0/H)). The red 
lines are the best fits to the data, taking into account the errors on the metallicities. The yellow/blue/green regions represent the Icr 
errors associated to the gradients. Blue label indicates the galaxies for which the gradient is positive within Icr, the green ones are 
those for which it is negative within the same limits. For each galaxy we have indicated the dynamical (Rot=rotating disk, NRot=no 
rotation) and environment (lso=isolated, lnt=interacting) classes. 



posed several scenarios to explain their discovery: (i) a radial 
redistribution of the metal-rich gas produced in the nucleus, (ii) 
supernovae blowing out metal-rich gas, enriching the IGM, then 
falling onto the outer parts of the disk, (iii) the result of a past 
interaction. 



At high redshift, Cresci et al.| (|20T0| have recently studied 
with SINFONI the metalUcity distribution of three Lyman-break 
galaxies at z ~ 3 in the AMAZE/LSD sample. They were able 
to derive metallicity maps for the three galaxies with high SNR. 
In each case, they discovered a positive gradient, comparable to 
the ones we have found. They favour the scenario in which these 
positive gradients would be produced by the infall of metal-poor 
gas into the center of the disks, diluting the gas and lowering its 



metallicity in the central regions. The authors claim that the dis- 
covery of positive gradients in high-redshift disks, pre-selected 
to be "isolated", is a direct evidence for cold gas accretion as a 
mechanism of mass assembly. Such a conclusion could be bal- 
anced arguing that a merger remnant can keep, during a transient 
phase with a typical timescale of ~ 0.5 Gyrs, an inversed metal- 
licity gradient (eg. Perez et al.|201 1 Torrey et al.|201 l| i. 



In contrast, our study appears to show that among the seven 
detected positive gradients, only two galaxies are isolated, the 
others showing signs of interaction. Cold gas accretion toward 
the center of disks might thus not be the only process able to 
lower the central metallicity. 
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Table 3. Integrated metallicity (Z = 12 + log(0/H)), metallicity 
gradient, stellar and dark matter halo masses of the MASSIV 
"first epoch" sample galaxies. 



Galaxy 
(1) 

020106882 
020116027 
020126402 
020147106 
020149061 
020164388 
020167131 
020182331 
020193070 
020208482 
020214655 
020217890 
020239133 
020240675 
020255799 
020261328 
020278667 
020283083 
020283830 
020294045 
020306817 
020363717 
020370467 
020386743 
020461235 
020461893 
020465775 
140083410 
140096645 
140123568 
140137235 
140217425 
140258511 
140262766 
140545062 
220014252 
220015726 
220071601 
220148046 
220376206 
220386469 
220397579 
220544103 
220544394 
220576226 
220578040 
220584167 
220596913 
910193711 
910279515 



(2) 
^.73 
^.42 

^.49 

^.58 

^.89 
^.76 



;.80 
;.39 

;.74 
;.6l 



;.57 
;.75 
;.53 
;.78 



8.72 
9.23 
8.75 
8.65 
9.02 
8.95 
8.69 
8.69 
8.78 
8.75 



!.40 
!.32 
!.38 
!.89 
!.55 
!.75 
!.64 
!.53 

!.77 
!.60 



A,.Z 
[dex kpc"'] 

(3)_ 

0.037 ± 0.026 
0.03 ± 0.024 



0.04 + 0.053 10.10 



0.0035 + 0.015 



0.045 ± 0.026 



0.0046 ± 0.039 
-0.022 ± 0.054 
0.0099 ± 0.024 
-0.078 ± 0.023 



-0.02 + 0.041 

0.019 + 0.074 
-0.051 +0.055 
0.023 + 0.0035 
-0.061 +0.037 
0.098 + 0.11 
-0.033 + 0.012 
0.021 + 0.0088 
0.00041 + 0.02 



0.04 + 0.021 

-0.022 + 0.02 
0.0063 + 0.012 

0.034 + 0.022 
-0.0073 + 0.01 
0.0031 +0.013 

-0.03 + 0.013 

-0.0091 +0.028 



log(M*) 
We] 
(4) 

Q QQ+().22 

10 09+" "' 

10.09+°! 

,+o:t3 



10.18! " 
10.13!p 
10.08!;' 
10.72+°" 
10 15+'^-^'^ 
10 17+""' 
10 02+°"' 

Q QQ+0 19 

9 89+° 
9 96+° 
9 87!° '^ 
10.01+^1 
10.28!?;« 
10 05+°^' 
10 37+° " 
9.80!" 
9 76+'^-^^ 
9 68+" " 
10.57!°;1^ 

9.88!SiS 
10 36+° 

9 66+8:^5 
-Oil 

10 12+"^" 

-S-2Q 
10.07!;! 

10.40!«;i 

q 7^+0.39 
^•'-^-0.39 

10.07!!J-2« 
10.84!";" 
10.80!«;i^ 
9 84+""' 
10.60!;! 
10.78!«;« 
10.77!°;^' 
10.81!°;« 

11 22+°" 
10.67!°;^^ 
10.80!S;1^ 
10 23+" '" 

10-71 

10 34+°-^3 
10 31+" 2' 

lU.Jl 023 

10.72+"" 
11.21 
10.68 

Q QQ+0.42 

10.79!":i 



log(Mh) 
[Me] 

(5) 

11.59+»f 
9 48+°-'^ 

Q 41+2 87 

-12,27 

11.33 



m 
i\ 

-0.42 

1 1 92+°-2' 
10.46+"? 



11.58 
11.57 
11.52 



11.84!°J^ 

10.26!38, 

11 41!" '2 

10.95+i^f 
10.56!°-^2 

12 02+° '^ 
-Q-2Q 
J 2 43+0.23 



10 12+^*'' 
10 35+^ 
10 15+°-2* 

ii.07!°:5| 

10.60!°! 
10 72+°-^' 

(J ^7+2 41 

12.75+t^^ 



!0 



10.43 j3 28 

10 67+° '' 

^"•"'-0.36 

12.85+°'6 

1 1 54+°-22 

1 1 42+2 86 
ii-t^_14 27 

I 9 9 Q +0-24 

-Q-5S 

I I 57+°-23 

12 33+^'*' 



10.38!^;^'^23 
12 17+0-'5 

^^•^ '-0.24 

1 1 30+°-^' 

I l.JU 54 

g Jg+I-02 

I I 62+''-^^ 
10 54+°-^' 

lU.Jt 044 

Q 77+0.49 
^- ' -10.09 

12.48!|,««3 
12.31+°^" 



-Ml 



11.69_ 
10 56+**-^' 

lU.JU 042 

12.48+!;?* 



3.2. Relations with global galaxy properties 

Several authors found correlations between gas-phase metallic- 
ity gradients and other global physical parameters of galaxies. 
The observed metallicity gradients are a function of the mor- 
phological type ( Vila-Costas & Edmunds||T992 Marquez et al. 
2002 1: they are steep in late- type spirals and almost flat for early- 
type spirals. Further, in the local Universe the absolute value of 




RIS NRIS RIN NRIN 



Fig. 6. Left: Metallicity gradients distribution for 26 MASSIV 
galaxies. The histogram in black line represents the whole sam- 
ple, the red and dashed blue ones respectively the isolated 
and interacting galaxies distribution. Right: Metallicity gradi- 
ents with respect to the kinematic al type (from left to right. 
Rotating Isolated, Non-Rotating Isolated, Rotating Interacting, 
Non-Rotating Interacting). 




40 60 
fjy [km.s-^] 
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Z - 12 + log(0/H) [dex] 



Fig. 7. Left: Metallicity gradient versus the mean velocity disper- 
sion of MASSIV galaxies. The black squares are the individual 
galaxies, and the red ones represent the median values for 3 bins 
of cr,,, along with the standard deviation in each bin represented 
by the error bars. Right: Metallicity gradients versus the inte- 
grated metallicity of each galaxy, the red line is the best fit to the 
data. 



the gradients seems to decrease with increasing luminosity (less 
luminous galaxies have steeper metallicity profiles) as predicted 
by mod eling (Prantzo s & Boissier 2000) and verified by obser- 
vations ( |Garnett et al.|I997t|van Zee et al.|1 998). 

Our high-redshift sample did not verify the latter correlations 
and we were not able to test a correlation with the morpholog- 
ical type, as we do not have that information for our MASSIV 
sample. The right panel of Fig. |6] shows no clear trend between 
the strength of the metallicity gradient and the kinematical type, 
nor any correlations with their close environment (isolated or 
interacting). We can however notice that i) among the isolated 
objects, the non-rotating galaxies have on average flatter gradi- 
ents than rotating disks, and ii) the fraction of positive gradients 
is higher in interacting systems compared with isolated galaxies. 

One of the two weak correlations that might be present in our 
sample is shown in the left panel of Fig. |7] The strength of the 
gradient seems to correlate with the velocity dispersion of the 
galaxy. The latter is derived on beam smearing corrected veloc- 
ity dispersion maps obtained after velocity field modeling (see 
Epinat et al.|20I I[ ) and thus reflects the true velocity dispersion 
of the gas. 

This apparent correlation seems to be driven by the fact that 
galaxies with high gas velocity dispersion show shallower, of- 
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Fig. 8. Mean infall rates over 4 Gyr versus the mass of the host- 
ing DM halo of each MASSIV galaxies. The objects classified as 
rotators are circled in black. On the lower left comer is displayed 
a typical errorbar on the DM halo mass. 



ten positive, metallicity gradients. This correlation has, to our 
knowledge, not been observed previously in the local universe, 
where the velocity dis persion in late-typ e objects is usually low 
with cr,, ~ 20 km s ' ( jEpinat et al. 2010 1. Could turbulent phys- 
ical conditions in the ISM of high-redshift galaxies be at the ori- 
gin of the shallow, sometimes positive gradients? This question 
is difficult to address here, considering the relatively low spatial 
resolution of our data and the scatter in Fig. |7] At face value, 
the positive gradients in our z ~ 1-2 galaxies might be related to 
the perturbed physical conditions/motions in the ISM of high-z 
galaxies, as opposed to the continuous metallicity gradients ob- 
served in the relatively quiet ISM of the local spirals. 

Finally, we observe as well a tentative anti-correlation be- 
tween the metallicity gradient of each galaxy and its integrated 
metallicity (see Fig. [7] right panel). Metallicity gradients are 
more frequently negative in metal-rich galaxies and more fre- 
quently positive in low-metallicity galaxies. If real, this be- 
haviour would support the scenario in which infall of metal- 
poor gas from the IGM into the center of the disks drives the 
positive gradients. This infall of pristine gas would be able to 
reverse the gradient by diluting the central gas metallicity, and 
lowering the overall metallicity of the galaxy at the same time. 
The main question would remain: where does the metal-poor gas 
come from? Accretion of cold gas from the DM reservoir and/or 
interaction-triggered gas infall/capture from companions? 



3.3. Gas infall rates 



As suggested already by previous studies (eg. Werk et al.|2010 



|Rupke et al.||2010) , during an interaction, metal-poor gas from 
the outskirts of the galaxy could radially flow towards the center 
and dilute the metallicity in the inner, high star-forming regions. 
At the same time, metal-enriched gas could be transported to the 
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Fig. 9. Mean infall rates versus the star formation rates derived 
with the Ha luminosity and corrected for extinction. The line is 
the best fit to the data. 



outer parts during the interaction, overall resulting in a flattening 
of the metallicity gradient. 

The two scenarios explaining the positive metallicity gradi- 
ents - cold gas accretion and gas redistribution during interac- 
tions - have in common that metal-poor gas needs to be trans- 
ported efficiently to the center of the objects on time scales 
shorter than the star formation. 

We consider the possibility of infall of pristine gas onto the 
disk in the context of a chemical evolution model in order to 
explain the positive gradients. Our toy model for the chemical 
evolution of galaxies assumes (i) an instantaneous recycling ap- 
proximation (IRA), (ii) the infall of metal-free gas onto the disks, 
and (iii) radial flows of gas into the disks. This leads to the fol- 
lowing equations: 



— + ^(rvZ) = y{\ - Rm) - z(l - R)ilf{t) 
at ror 



(3) 
(4) 



where g is the gas surface density, ij/ is the star formation 
rate per unit area, R is the "returned fraction", / is the infall rate 
(per unit area) of intergalactic metal-free gas, Z = is the metal 
content, y is the stellar yield, and v is the velocity of the radial 
gas flow (positive toward the outer parts). Following a Schmidt- 
Kennicutt law for the star formation rate density (E^f, = Cg", 
C - 1.6 X 10"^^, n - 0.7) we solve these equations for each 
galaxies in our sample, assuming typical values for C, n, y, v and 
R (see above and below). After combining equations ([3]) and Q, 
we end up with the following equation allowing to derive the 
infall parameter /: 



dz dz 



■■y(l-R)C''"^ 



\ln.\-\in , 



■zC^'"<p-^'"f 



(5) 
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Assuming the physical parameters involved in eq. (|5]l to re- 
main constant for over ~ 4 Gyr (i.e. valid for galaxies at z ^ 1 .2), 
we derive a mean infall rate per unit area that can be integrated 
since the formation epoch. To this end, we use a true yield value 
of y - 0. 019, and a returned fraction of ~ 40%: R - 0.4 
(Erb 2008[ l. The radial flow velocity of high-redshift galaxies 



Complementary analyses focused on the kinematical classifica- 
tion and the fundamental scaling relations can be found in two 
companion papers ( Epinat et al.||201 1\ Vergani et al.||20lT i 



is known to be higher than in local spirals (e.g. ~ 1 km s in 
the case of the Milky Way). This velocity can be approximated 
to be the velocity towards the center of the migrating clumps - 
- in which the star formation ta kes place - which then leads to 
velocities of 10 - 50 km s ' (Bournaud et al. 2007| l. For our pur- 
pose, we used for each MASSIV galaxy the estimated velocity 
dispersion reported in Epinat et al.,(2011j ). 

For the star formation rates, we used those derived from the 
Ha luminosity, corrected for extinction, and listed in Table [2] 
Although the Ha luminosity only reflects the instantaneous star 
formation rate, we assume it to be an average value since the for- 
mation epoch. As our galaxies were selected to be star-forming, 
this is most likely an upper limit on the average star formation. 

The so-calculated global infall rates range from a few to sev- 
eral hundreds Moyr These values are high compared to the 
output of hydro-dynamical simulations of typical disk galaxies 
(50 Moyr"'). However, the rate of infalling gas is believed to de- 
pend on the size and mass of the host DM halo, being able to 
reach a few hundreds of Mgyr"' for the most massive galaxies 
( |Keres et al.|[2005j |Erb„2008) . We have therefore derived DM 
halo masses for each galaxy in our sample (see Section [23] l and 
compared it to our calculated infall rates (see Fig. [8]). As ex- 
pected there is a tendency for galaxies to accrete more gas in the 
most massive halos. We further divided the sample in three bins 
with respect to the (extinction-corrected) star formation rate. In 
that case, we noticed that galaxies with the highest SFR also have 
the highest gas-infall rate, whereas galaxies with the lowest SFR 
show infall rates typical for local galaxies (~50 Mgyr"'). 

Two MASSIV galaxies deviate from the global trend in 
Fig. [8] VVDS220397579, located at the left end of the graph, 
has a very low DM halo mass but an important infall rate 
(~ 650 Moyr-i), while VVDS220376206 has by far the highest 
infall (> lO"* Moyr"') and star formation (> 250 Moyr"') rates . 
For VVDS220397579 it is very likely that its DM halo mass has 
been underestimated. Indeed, this galaxy could be a face-on disk 
as the kinematical modeling returned a very low rotation veloc- 
ity of ~ 9 km s"'. As the DM halo mass is directly proportional 
to the maximum rotation velocity (see equation [TJ the DM halo 
mass for VVDS220397579 is most probably only a lower limit. 

In Fig.|8] we have circled the objects which, according to our 
kinematical classification, are isolated and rotating. For those 
objects the DM halo mass is more reliable, and they appear in- 
deed to have the largest DM halo masses. In order to investigate 
the relationship between the infall rate and the star formation 
rate, we show in Fig. [9] the mean infall rates against the SFRs 
corrected for extinction. As expected from equation |5] a clear 
linear correlation between these two quantities can be seen: the 
infall rates appear to be directly proportional to the star forma- 
tion rates. The best fit to the data gives a slope of 2.71 for this 
relation, which is defined as the infall parameter fi. This was al- 
ready shown to be a typical value for galaxies at these redshifts 
( |Erb|2008| l. 



4. Summary and conclusions 

This paper presents the chemical abundance analysis of the 
"first epoch" sample (50 galaxies) of the MASSIV survey. 



well as a description of the whole MASSIV sample ^Contini| 
|etal.|2011D . 

We have been able to measure emission lines ([Nii]6584 
and Ha) in 34 integrated spectra of MASSIV galaxies. Within 
this sample, we have identified one galaxy hosting an AGN. 
Chemical abundance estimates could be obtained for the remain- 
ing 33 star-forming galaxies. 

For 26 galaxies, we have been able to derive a metallicity 
gradient, defining annular regions around the peak of the Ha 
flux (which corresponds, in most cases, to the kinematical cen- 
ter). While just over half of our sample is compatible with a zero 
metallicity gradient, about a quarter of the galaxies shows pos- 
itive gradients: the metallicity increases from the center to the 
outer parts of the galaxy. Among these latter (seven) galaxies, 
four are classified as interacting systems, one is probably a chain 
galaxy, and two are classified as isolated. 

Flat, more rarely positive, metallicity gradients have already 
been found in interacting galaxies in the local Universe. They are 
explained by the infall of metal-poor gas onto the central parts 
during the encounter ( Werk et al.|2010[|Rupke et al.|2010l l. 

Two of our galaxies showing a positive metallicity gradient 
are classified as isolated and do not show any sign of recent in- 
teraction. Three analogue objects were reported as isolated disk 
galaxies at z ~ 3 by Cresci et al. ( 2010| l. In these cases, cold 
gas accretion onto the central regions of the disks seems to be 
the most plausible scenario. Even if tentative, there is a very 
interesting trend as a function of redshift. At z ~ 3 almost all 
isolated galaxies have a positive gradient, whereas this fraction 
drops to ~ 15 - 20% at ~ 1.2 and is almost equal to zero in the 
local universe. If cold accretion is the main process to explain 
the positive metallicity gradients in isolated disks, that suggests 
that the epoch where cold accretion dominates the mass assem- 
bly processes is at z > 2. 

We noticed in our sample a tendency for galaxies with the 
highest gaseous velocity dispersion to have a shallow/positive 
gradient, which highlights the different physical conditions ob- 
served in the ISM of high-redshift galaxies (high velocity dis- 
persion compared to local spirals). We discovered a weak cor- 
relation between gradient and global metallicity of the galaxy: 
metal-poor galaxies preferentially have flat/positive gradients 
while metal-rich ones tend to display negative metallicity gra- 
dients. This behaviour can also be explained by the infall of 
metal-poor gas onto the center of the disks, diluting the overall 
metallicity. 

Finally, applying a simple chemical evolution model with ra- 
dial flows of gas, we estimated infall rates of pristine gas onto the 
disks. We found values up to several hundred of Mo per years, 
and a tendency for the maximum infall rate to increases with the 
DM halo mass. 

The analysis of the spatially-resolved metallicity of galaxies 
will be further extended to the full MASSIV sample, enabling 
a better statistics and hence a stronger interpretation in terms of 
galaxy assembly scenarios. The final sample will also allow to 
establish the mass-metallicity and fundamental metallicity rela- 
tions at 0.9 < z < 1.8, adding constraints on the evolutionary 
status of star-forming galaxies at high redshifts. 
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